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FOREWORD 


This  report  presents  a  functional  and  mathematical  description  of  the  ground 
station  alert  predictor.  The  alert  predictor  generates  the  times  during  which  a 
satellite  will  be  visible  to  a  ground  tracking  station  and  provides  ground  antenna 
pointing  information.  This  is  an  invaluable  tool  that  is  useful  for  scheduling 
ground  station  operational  activities.  The  methods  used  to  generate  this  informa¬ 
tion  are  computationally  quite  efficient  so  that  long-term  predictions  can  be  done 
rapidly.  Several  new  innovations  are  included  in  the  algorithms:  the  effects  of 
atmospheric  drag  upon  satellites  are  taken  into  account,  and  a  satellite  orbit 
adjust  may  be  included  during  the  prediction  interval. 

Released  by: 


0.  F.  BRAXTON,  Head 
Strategic  Systems  Department 
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SECTION  1 
INTRODUCTION 

The  primary  function  of  the  ALERT  package  is  to  predict  the  times  during 
which  a  satellite  will  be  visible  at  a  ground  tracking  station  and  provide 
approximate  ground  antenna  pointing  information.  The  methods  used  to  gen¬ 
erate  this  information  are  computationally  quite  efficient,  so  that  the 
in-view  times  and  associated  information  for  a  given  satellite  and  ground 
station  can  be  rapidly  predicted  for  fairly  long  prediction  intervals 
(  ~  30  days).  An  additional  optional  feature  exists  which  accounts  for 
the  inclusion  of  a  satellite  orbit  adjust  during  the  prediction  interval. 

A  functional  overview  of  the  ALERT  generator  is  shown  in  Figure  1-1.  As 
can  be  seen  from  this  figure,  the  ALERT  generator  is  comprised  of  four 
basic  computational  functions:  the  process  flow  supervisor  (PFS),  the 
Brouwer  mean  element  converter  (BMC),  the  satellite-station  in-view  pre¬ 
dictor  (SIP),  and  the  Brouwer  orbit  generator  (BOG).  Also  shown  on  this 
figure  is  the  inter-function  data  flow.  A  detailed  description  of  each  of 
these  functions  is  presented  in  the  following  sections. 

As  user  supplied  input,  the  ALERT  generator  requires  the  station  and  satel¬ 
lite  and  associated  frequency  information  for  which  in-view  periods  are  to 
be  predicted,  as  well  as  the  start  (tg)  and  stop  (tg)  times  of  the  time  interval 
over  which  predictions  are  to  be  made.  Also  included  in  the  input  is  the 
option  selected  by  the  user  for  the  inclusion  of  the  impact  of  orbit  adjusts 
occurring  during  the  prediction  interval  on  in-view  periods.  In  order  to 
initiate  its  processing  the  ALERT  generator  requires  the  following  additional 
data  as  shown  on  Figure  1-1:  station  latitude,  longitude,  and  height;  mean 
Brouwer  elements  and  decay  rates;  and  post  orbit  adjust  epoch,  inertial 
position  and  velocity  vectors  if  the  effect  of  an  orbit  adjust  is  to  be 
included  in  the  prediction  span.  The  predicted  inview  periods  and  associated 
station  antenna  pointing  geometry  that  were  computed  for  the  prediction  interval 
are  written  to  hardcopy  at  the  end  of  the  computational  cycle. 


POST  ORBIT 
ADJUST  DATA 


*OA-  '<t'OA>-  r(t'OA) 


PREDICTED  IN-VIEW  PERIODS  AND 


^ANTENNA  POINTING  GEOMETRY 


V> 


-y  *  •  . 

' 0  a 


FIGURE  1-1.  ALERT  FUNCTIONAL  OVERVIEW  END  DATA  FLOW 


SECTION  2 


THE  PROCESS  FLOW  SUPERVISOR  (PFS) 

2.1  FUNCTIONAL  DESCRIPTION 

The  principal  tasks  performed  by  the  PFS  are  to  receive  and  retrieve  data, 
direct  processing  flow,  and  output  computed  results.  Specifically  the  PFS 
function: 

(i)  receives  input  data; 

(ii)  retrieves  station  coordinates; 

(iii)  retrieves  initial  SPASUR  Brouwer  mean  elements  and  decay  rates; 

(iv)  receives  predicted  orbit  adjust  times  and  post  adjust  vectors; 

(v)  cycles  computations  through  BMC  if  required; 

(vi)  cycles  computations  through  SIP;  and 

(vii)  receives  SIP  predictions  and  outputs  them. 

The  flow  of  the  PFS  function  is  presented  in  Figure  2-1. 

It  should  be  noted  that  tasks  (iv)  and  (v)  above  are  optional.  Should  the 
user  elect  to  use  one  of  the  two  orbit  adjust  options  and  an  orbit  adjust 
falls  within  the  prediction  span,  the  PFS  retrieves  the  predicted  orbit 
adjust  time  (tQA)  and  the  associated  predicted  post  adjust  position  and 
velocity  at  a  time  t^A  very  near  (but  later  than)  tQA,  r(tQA)  and  r(tQA).  These 
osculating  Cartesian  vectors  are  passed  to  the  BMC  function  where  they  are 
converted  to  Brouwer  mean  elements  for  the  time  tQA. 

As  noted  above,  the  orbit  adjust  option  can  operate  in  two  modes.  If  the 
orbit  adjust  option  flag  IOAOP  is  set  to  zero,  the  option  is  not  used  and 
SIP  predictions  are  generated  using  only  the  initial  element  set  obtained 
from  the  SPASUR  data.  However,  if  it  is  set  to  one  or  two,  the  retrieval 
and  conversion  of  the  post  adjust  vectors  take  place,  as  described  above. 

When  IOAOP  is  set  to  one,  the  prediction  span  is  divided  into  two 

partitions  using  tgA  as  the  boundary,  and  the  initial  Brouwer  elements  retrieved 
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FIGURE  2-1.  PROCESS  FLOW  SUPERVISOR  LOGIC  FLOW 
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FIGURE  2-1.  PROCESS  FLOW  SUPERVISOR  LOGIC  FLOW  (CONTINUED) 
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2.1  FUNCTIONAL  DESCRIPTION  (con't) 


by  the  PFS  from  the  SPASUR  element  data  are  used  for  SIP  in-view  pre¬ 


dictions  until  tg^  is  reached. 


For  times  later  than  tg^  the  post-orbit 


adjust  Brouwer  mean  elements  generated  by  the  BMC  function  are  used  for 


SIP  predictions.  When  IOAOP  is  set  to  two,  SIP  predictions  are  generated 
only  for  times  later  than  tg^  using  the  Brouwer  element  set  generated  by 
the  BMC  function. 


2.2  SPECIAL  FEATURES 

Several  points  associated  with  the  orbit  adjust  option  are  worthy  of  note. 
First,  the  orbit  adjust  options  (IOAOP  t  0)  may  be  selected  only  for  those 
satellites  which  are  supported  by  the  PULSAR  system1.  If  the  option  has 
been  selected  for  a  satellite  other  than  those,  the  option  is  denied  and 
processing  continues  as  if  the  option  had  never  been  selected.  Secondly, 
if  the  orbit  adjust  option  has  been  validly  selected,  the  PFS  uses  the 
decay  rates  obtained  from  the  SPASUR  element  file  for  the  epoch  tg  as  the 
dec  .y  rates  for  times  later  than  tg^. 

The  epoch  tg  above  corresponds  to  the  time  of  the  Brouwer  mean  element  set 
obtained  from  the  SPASUR  data  nearest,  but  not  later  thar,  the  start  time 
of  the  prediction  interval  tg,  i.e.  tg  <_  tg.  If  an  orbit  adjust  occurs 
within  the  time  interval  between  tg  and  tg  processing  is  still  ’-erformed 
as  described  above.  However,  the  alert  schedule  only  for  the  time  inte-val 

f  tB’  tE  1  be  outPut- 


Parks,  A.O.  and  Hicks,  T.I.,  A  Mathematical  Description  of  the  PULSAR 
Doppler  Tracking  Data  Editor,  NSWC/DL  TR  82-591.  Dahlgren, 


Virginia,  1982. 
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2.3  PROCESSING  EQUATIONS 


The  semi -major  axis  read  from  the  SPASUR  data  is  the  Kaula  semi -major 
axis  ak  expressed  in  earth  radii.  This  is  converted  in  the  PFS  to  the 
Brouwer  mean  semi -major  axis  a"  via  the  transformation 


a.  a 
k  e 


(2.1) 


where 


3J,(1  -  3/2  sin  -4" ) 


„2v3/2 


4  a^  (1  -  e"2) 


In  the  above  expressions  ag  is  the  earth's  semi -major  axis;  J2  is  a  zonal 
harmonic  gravitational  constant;  and  l "  and  e"  are  the  Brouwer  mean  inclina¬ 
tion  and  eccentricity,  respectively. 

Similarly,  the  semi-major  axis  decay  rate  obtained  from  the  SPASUR  data  is 
the  time  rate  of  change  of  the  Kaula  semi -major  axis  ak-  This  is  converted 
in  PFS  to  the  time  rate  of  change  of  the  Brouwer  mean  semi -major  axis  a" 
through  application  of  the  expression 


5"  =  ^(THr)273  +2akae*  [d*2X)(l 


(2.2) 


(2.3) 


where 


■!■(*-)•■  -  4)1  - 


In  the  last  expression  it  has  been  assumed  that 


(2.4) 


(2.5) 


Values  for  e"  and  the  rate  of  change  of  the  mean  motion  n  are  also  obtained 
from  the  SPASUR  data. 
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SECTION  3 

THE  STATION-SATELLITE  INVIEW  PREDICTOR  (SIP) 

3.1  FUNCTIONAL  DESCRIPTION2 

The  station-satellite  in-view  times  can  be  obtained  by  considering  the 
vertical  distance  A  of  the  satellite  above  a  horizontal  plane  tangent  to 
the  earth  at  the  station's  location.  An  inview  period  is  considered  to 
begin  or  end  when  the  satellite  is  in  the  horizontal  station  plane,  i.e., 
when  A  =  0.  Since  the  earth  rotates,  the  satellite  will  dip  below  the 
zero  level  and  may  not  appear  above  the  horizontal  plane  for  many  revolu¬ 
tions. 

To  initiate  SIP  processing,  it  must  first  be  determined  where  the  satellite 
is  at  time  tg  with  respect  to  the  maximum  points  of  the  function  describing 
A.  Once  this  is  done,  the  first  maximum  is  approached  in  fractional  parts 
of  a  satellite  period  P  until  a  close  proximity  of  the  maximum  A  is  reached. 
This  is  accomplished  by  using  the  signs  and  magnitudes  of  the  A  and  its 
first  and  second  time  derivatives.  A1  and  A",  to  determine  how  far  to  step 
forward  in  time.  When  the  maximum  A  has  been  approached  sufficiently  close, 
and  it  is  greater  than  zero,  the  exact  time  of  maximum  A  is  obtained  itera¬ 
tively.  This  time  is  used  to  compute  approximate  rise  and  set  times,  which 
are  then  refined  using  an  iterative  method.  This  procedure  is  repeated  to 
find  successive  maximum  As  and  associated  rise  and  set  times  until  the  end 
of  the  prediction  interval  t£  is  exceeded. 

For  each  in-view  period  geometric  parameters  are  computed  and  stored.  When 
the  end  of  the  prediction  interval  is  exceeded,  the  in-view  times  and  associ¬ 
ated  geometry  computations  are  passed  to  the  PFS  for  output.  The  process 
flow  of  SIP  is  presented  in  Figure  3-1. 


2Castro,  H.E.,  A  New  Method  for  Finding  Station  Prediction  Alerts,  NSWC/DL 
TM  K-l/63,  Dahlgren,  Virginia,  1963. 
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FIGURE  3-1.  STATION-SATELLITE  IN-VIEW  PREDICTOR  PROCESS  FLOW 
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3.2  PROCESSING  EQUATIONS 

The  orbital  period  P  is  computed  using 


n  2tt  3/ 

P  =  - a  2 


(3.1) 


where  y  is  the  gravitational  constant  and  a  the  semi -major  axis. 

The  vertical  distance  A  of  a  satellite  from  a  horizontal  plane  tangent  at 
a  station  location  is  given  by: 


A  =  cos<|>  (x  cosa  +  y  sin  a)  +  Zsint  -  RQ  , 


(3.2) 


where 


w(t  -  tye)  +  X  , 


Ro  =  ag  (l-e2  sin2#)*4  +  h  , 

<j>  =  station  geodetic  latitude, 

X  =  station  longitude, 

w  =  earth's  rotation  rate, 

a  =  earth's  equatorial  radius, 

e 

tve  =  universal  time  of  transit  of  vernal  equinox, 
e  =  eccentricity  of  earth, 

x,  y,  z  =  satellite  inertial  position  components,  and 
k  =  height  of  station  above  reference  ellipsoid. 
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3.2  PROCESSING  EQUATIONS  (con't) 


The  first  and  second  time  derivatives  of  A  are  given  by 


A'  =  cos  4>  |  [y-u>x)  sina  +  (x  +  wt/)  cosa  J  +  z  sin<f>  , 


(3.3) 


A"  =  cos  <J>  [({/  -  2a)X  +  sina  +  (x  +  2wy  -  w2x)  cosa  ]  +  zsin<j>  ,  (3.4) 


where 


x,  y,  z  =  satellite  inertial  velocity  components. 


x  =  -  y  — —  , 
r3 


t/  =  -  y  y  *  and 


z  =  -  y 


The  expression  required  for  the  iterative  solution  for  the  time  of  maximum 
A,  i.e.  tm,  is  given  by 


tUn>  .  ttt)  .  ‘L^l  , 

m  m  A"(t(/L}) 

m 


(3.5) 


where  the  superscripts  in  parentheses  indicate  ^  s  obtained  from  the  x 
iteration.  Equation  (3.5)  is  assumed  to  have  converged  when  the  tolerance 
condition 


tU+l)  _  U)\  <  T 

m  m  1  i  * 


(3.6) 


W.v 


-fcwTV  W  ITT  rr  VHJ  XT>  UT  V  -  I  i.  VTLiT*  iT*U^  *-fc-  »  im  i.  >  l  lit 
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3.2  PROCESSING  EQUATIONS  (con't) 


is  satisfied.  The  first  estimate  of  tm  to  be  used  in  equation  (3.5)  for 
tj°)  is  the  value  of  t  obtained  prior  to  entering  the  A  section  of  logic 
in  Figure  3-1. 

First  estimates  of  the  rise  time  tr  and  set  time  t$  associated  with  the  time 


of  maximum  A,  t^,  are  given  by 


K  = 

r  m 


JUT 

"  1  A" 


and 


tc  =  tm  + 
s  m 


V? 


(3.7) 


These  estimates  are  used  to  initiate  iterative  solutions  for  more  exact 
rise  and  set  times  using  equation  (3.5).  Convergence  is  assumed  to  have 
occurred  when  equation  (3.6)  is  satisfied  for  both  tp  and  ts. 

3.3  1N-V1EM  GEOMETRY  PARAMETER  COMPUTATIONS 

The  time  of  maximum  A,  i.e.,tm,  is  used  to  initiate  the  iterative  solution 
for  the  time  of  closest  approach  t^,  via  the  following  expression: 


.(x>l)  _  .(>c) 

lCA  TA 


A  - 


A2 


(3.8) 


r  -  r 


t  =  t^ 
CA 


where 


A  =  (x  -  a,  cosa)  (x  +  fl,  ui  sina)  +  (y  -  olx  sina)  {y  -  ax  u>  cosa)  (3.9) 


+  (z-b^z  , 


A  =  (x  -  at  cosa)  (x  +  a ^  u>2  cosa)  +  ix  *  ai  wsina)2  +  (y  -  sina) 


••  #  ••  # 

( y  +  a i  w2  sina)  +  {y  -  ai  wcosa)2  +  (z  -  bj  z  +  z2  , 


(3.10) 


|r  -  rQ|2  =  (x  -  aj.  cosa)2  +  (y  -  ax  sina)2  +  (z  -  bx)2  , 

-12- 


(3.11) 
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3.2  PROCESSING  EQUATIONS  (con't) 


and 


ai  =  |ae^  "  £2  s^n2<^)"Js  +  cos<f>  » 


b1  =  ag(  1  -  e2)  (1  -  e2  sin2<j>)~a  +  h  sin<f> 


The  azimuths  A  at  tr,  t CA>  and  t  are  computed  using  the  equation: 


(3.12) 


(3.13) 


A  =  tan-1 


-(x  sin  a  -  ycos  a) 


-sin<p  (xcosa+  t/sina)  +  cos<j>  |  (z  +  ag  e2  sin<j>)/(l-e2  sin2<j>)SsJ 


(3.14) 


If 


x  sin  a  -  y  cos  a  |  <  t 


(3.15) 


and 


|sin<|>  (xcosa  +  y  sin  a)  -  cos  cp  J(z  +  age2  sin<|>)/(l  -  e2  sin2^)*5  J  |  <  (3.16) 


then  A  is  set  to  zero. 


The  Doppler  frequency  An  (scaled  to  100  me  and  including  the  offset  fre¬ 
quency)  is  computed  at  tf  and  ts  using  the  relation: 


An 


=  ( - - - ^  f  (x-a,  cos  a)(x  +  a,  w  sina)  -  ( y-a ,  sina)(y  -  a.  wcosa) 

V  .2997928  /  l  1  11 

+  (z  -  bjz  ]  •  [  (x  -  cosa)2  +  (i f-a1  sina)2 


(z-  b  )2  j  **(  —  )+  lOOf 

i  1  V  108  /  os 


(3.17) 
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3.2  PROCESSING  EQUATIONS  (con't) 


where  f  is  the  satellite  frequency  and  fQS  is  the  satellite  off-set 
frequency. 


The  elevation  E  at  t ^ 


is  computed  using  the  expression 


and  the  slant  range  p  at  tCA  is  obtained  from 

[  (x  -  cosa)2  +  (y  -  sina)2  +  (2  -  bj)2]*5 


P  ■ 
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SECTION  4 

THE  BROUWER  ORBIT  GENERATOR  (BOG) 

4.1  FUNCTIONAL  DESCRIPTION 

The  Brouwer  orbit  generator  function  is  used  by  the  station  alert  pre¬ 
diction  software  package  to  supply  predicted  inertial  position  and  velocity 
vectors  for  the  satellite  being  processed.  These  vectors  are  used  by  the 
SIP  function  to  generate  station-satellite  in-view  periods.  The  bog  function 
is  also  used  by  the  BMC  function  to  convert  osculating  inertial  Cartesian 
coordinates  to  mean  Brouwer  elements  (it  is  assumed  in  this  case  that  the 
satellites  for  which  this  is  done  have  inclinations  which  are  not  in  the' 
neighborhood  of  the  Brouwer-Lyddane  singularity  at  -63.5°). 

Two  methods  are  available  for  generating  osculating  elements  in  the  BOG 
function.  As  long  as  the  satellite  inclination  is  not  in  the  neighborhood 
of  63. 5°,  the  Brouwer-Lyddane3*4  theory  modified  to  include  drag  effects  is 
used  to  generate  predicted  osculating  elements.  If  the  inclination  is 
sufficiently  close  to  63.5’,  a  Keplerian  update  method  using  secular  terms 
only  and  no  drag  effects  is  applied  to  the  mean  elements.  These  updated  ele¬ 
ments  are  then  used  as  the  osculating  elements  in  subsequent  SIP  calculations. 
In  each  case  when  BOG  is  being  utilized  by  the  SIP  function,  the  osculating 
elements  are  transformed  to  osculating  inertial  Cartesian  coordinates  before 
returning  to  SIP.  If  BOG  is  being  used  by  the  BMC  functions,  no  such  trans¬ 
formation  takes  place.  The  BOG  processing  flow  is  shown  in  Figure  4-1. 


3  Brouwer,  D.,  "Solution  of  the  Problem  of  Artificial  Satellite  Theory 
Without  Drag",  The  Astronomical  Journal ,  Vol.  64,  No.  1274,  1959, 
pp.  378  -  397. 

4  Lyddane,  R.  H.,  "Small  Eccentricities  or  Inclinations  in  the  Brouwer 
Theory  of  Artificial  Satellites",  The  Astronomical  Journal,  Vol.  68, 
No.  8,  1963,  pp.  555  -  558. 
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FIGURE  4-1.  BROUWER  ORBIT  GENERATOR  PROCESS  FLOW 
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4.2  PROCESSING  EQUATIONS  FOR  THE  BROUWER-LYDDANE  METHOD  (con't) 


r  =  n0t  J  1  +  y  Y2'n(3e2  .  1)  +  ^-  Y2'2n  [-15  +  16n  +  25n2 

+  (30  -  96n  -  90n2)e2  +  (105  +  I44n  +  25n2)e4| 

15  t  i 

+  —  Y/tie"2  [  3  -  3092  +  3504  ]  +  lQ"  +  nt2 


(4.2) 


3  3 

s"  =  n<)t  7  Y>'a  +  IF  V2  l'35  +  24n  +  25nZ 

+  (90  -  I92n  -  126n2)e2  +  (385  +  360n  +  45n2)04] 

+  [21  -  9n2  +  (-270  +  126n2)e2  +  (385  -  189n2)04]+  9”  ,  (4.3) 

h"  =  n0t  -  3y2'0  +|y2'2  [(-5  +  12n  +  9n2)e  +  (-35  -  36n  -  5n2)e3) 


+  —  y4'  (5  -  3n2)e  (3  -  7e2)  +  h0" 


(4.4) 


The  long  period  (dependent  upon  9 ")  terms  are  computed  from: 


616  =  %"  V7  eH2r,2x  s1n*"sinsg"-I-LV  (3y2,23-  10Y4’Y)sin2g" 


35  _YV 
128  y* 


n2  X  sin,t"  sing"  +  —  ^V|y3'  +  — y5'  (4+3e"2)6 


T  [ 3y|-  e -  10Y4' y| 


siru."  sing"  + 


(4.5) 
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4.2  PROCESSING  EQUATIONS  FOR  THE  BROUWER-LYDDANE  METHOD  (con't) 


V  +  9'  =  9“  +  l"  +  % 


24Y, 


-  [-  3y2 1 2  I  2  +  e"2  -  11(2  +  3e"2)e2 


40(2  +  5eM2)e,,o"1  -  400e"2  06 


a'2! 


+  10Y4' 


Y,' 


2  +  e"2  -  3(2  +  3e"2)02  -  8(2  +  Se"2)©4^1  -8Oe"206a 

35  Y. 


r  V2  5 

4  s  'TT>Y 

sin  2g"  + 

5  r,3pn 


384  y2' 


n  e"  Xsin-c" 


35  Yr  ' 


--  X 


1152  Y2 1 


+  2e,l302  si  at" 


Y3e"02  5 

+ 


!I3q2 


-  e"(3  +  2e"2)  si  at"  + 


e"*0 


si  at" 


5  +  3202  a*1  +  8O04  a“2 


Y5 


4Y2'  siat"  64  Y2 


cos  3g" 

7  [-«"  -r~  (4  +  3e"2)  +  e"siat" 
L  siat" 


Y  * 

(26  +  9e"2)]  6  -  —  e"02  siat"  (4  +  3e"2) 


32  Y2. 


(3  +  1602a-1  +  4O04a"‘:)  +  h  — -  siru". 


Y 


( 


1+n3 


3  -  e"2 (3  -e"2 


>]♦£ 


e"(-32  +  81e"4) 

siat" 

4  +  3e"2  +  n(4  +  9e"2) 

5  V  2* 

—  n  6 
64  y2' 


cosg 


(4.6) 


TR  82-387 


4.2  PROCESSING  EQUATIONS  FOR  THE  BROUWER-LYDDANE  METHOD  (con't) 


35 Y5‘  e,,3e 

144  y2‘ 


j  nsin"1^1  +  sin-t"  [5  +  3202a_1  +  80e4cx“2 ] 


12  Y,' 


sin2g"  cosg"  +  r  j  -  3Y2'2  [ll  +  8O02a_1  +  2OO04cf 2 ] 

+  10Y,/  |  3  +  1602a  1  +  4O04a"2  j  sing"  cosg" 

35Y  1  . 

+ - —  e"30  Uxsin-U"  +  sini"  (5  +  320  V1  +  8O04cf2) 

576Y2 '  1  ' 

+  - — — -  Y  -  +  y  1  (4  +  3e"2)6  +  ~  Y  '  (4  +  3e"2) 

4y2 ' siac"  3  16  5  8s 


(3  +  1602a"1  +  4O04a-t)  sin2^."j  cosg"  . 


The  short  periodics  (dependent  upon  E',  ft  ,  £")  are  computed  from: 


a  =  a"t  +  a"  -  a"  — | -  (302  -  1)  + 


(1  -  e"  cos  E')3 


(4.7) 


[  302  -  1  +  3sin2-t"  cos(2g"  +  2^ 1 )  J 

e  =  e"  +  e"t  +  6  e  +  ~  —  — - —  j  3  -  e"2(3  -  e"2)^ 

2  n6  [l+n3  I 

+  3  +  e"  cosiT  •  (3  +  e"  coS(('' )  cosjf  +  3^--  " 

1  n6 

[  e"  +  j 3  +  e"  cosjj'  (3  +  e"  cos^')  cos<$ 1  j  cos  (26'  +  2g") 


(4.8) 
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4.2  PROCESSING  EQUATIONS  FOR  THE  BROUWER-LYDDANE  METHOD  (con't) 


2V  • 


iTY 


(1  -  02)  |  3  COS  (2gu  +  *•)  +  cos  (2gu  +  3<o] 


(4.9) 


-t  *  -t 


g-.9—  6  e  +  e"Y  '  0  si  rtf."  sin*'  s1n(2rf*  +  2g") 


2-4-  /«  l 


rrslitf' 


+  2e"Y2' eslru"  •  costf'  cos(2^'  +  29")  +  — y2'6  slrtf"  cos(2^'  +  2g")  (4.10) 


g  +  l  «  g*  +t'+-| 


-  6a(*'  -  l*  +  e"  sln^ ' )  +  (3  -  502) 

[  3  sin  (2<J'  +  2g")  +  3e"  sin  (2g"  +  f)  +  e"  sin  (2g"  +  3rf')]  J 


e"n2Y2 1 


4(1  +  n) 


2(302  -  1)  (a  +  1)  sin*1  +  3(1  -  02) 


[  (1  -  a)  sin  (2g"  +  *’)  +  (0+  1/3)  sin  (2g"  +  3*')]|  (4.11) 
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4.2  PROCESSING  EQUATIONS  FOR  THE  BROUWER-LYDDANE  METHOD  (con't) 


h  =  h'  +  [ 2e"Y2 *0  cos^'  +  y  y2'0  ]  sin(2g"  +  2^') 


-  e"  Y2'6  sin^'  cos(2|5'  +  29"}  -  3Y  '6  ({•  -  l "  +  e"sintf‘) 


and 


eM  =  — 


1  e"n3 


2  Y  ' 


ty*‘b-tVy 


sin  29" 


5  Yr* 


1  Y  ' 

—  — ~n3  sin*."  +  -^n3  sin*"(4  +  9e2")s 

4  Y  '  64  y2‘ 


35  Y  ' 

cos  9"  +  ToT"  n3e2"X  sin*"  cos  3g" 
384  Y  2 1 


-  —  y'n3  2(302  -  1)  (a  +  1)  sinrf*  +  3(1  -  e2)  (1  -  o] 


4.  ‘2 


+  (a  +  sin  (2 9"  +  3n  , 


(4.12) 


sin(2g"  +  6') 

(4.13) 
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4.2  PROCESSING  EQUATIONS  FOR  THE  BROUWER-LYDDANE  METHOD  (con't) 


the  eccentric  anomaly  E'  is  obtained  from  a  Newton-Raphson  iteration 
upon  the  Kepler  equation 


E'  -  e"  sin  E*  =  l"  , 


and  the  true  anomaly  is  found  from 


sin  |£'  = 


n  sin  E‘ 

1  -  e"  cosE' 


cos  i'  = 


cosE'  -  e" 


1  -  e"  cos  E* 


The  final  osculating  values  for  a,  and  k  are  computed  from  equations 
(4.8),  (4.10),  and  (4.12),  respectively.  Equations  (4.2),  (4.9),  (4.11) 
and  (4.13)  are  used  to  calculate  final  osculating  values  for  £,  g,  and  e 
for  the  following  relations: 


A  =  ecosf"  -  e6£sin£"  , 


B  =  esin£"  +  e&lcosl"  , 


l  =  tan"  (B/A) 


9  *  U  +  9)  -  l 


e  =  (A2  +  B2)*5 
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4.3  PROCESSING  EQUATIONS  FOR  THE  KEPLERIAN  UPDATE  METHOD 
The  Keplerian  update  method  treats  a,  e,  and  l  as  constants  and  updates 
lt  g,  and  h  using  only  their  secular  variations.  The  updated  elements 
are  given  by: 


and 


(4.22) 


n0t  [l  +  y  Vn  (3cosU"  -1)]  +  iQ"  , 


ynoY2't  (5  cos  2i"  -  1)  +  gQ"  , 


(4.23) 

(4.24) 


-  3n0  Y2't  cos  <c"  +  hQ* 


(4.25) 


The  constants  appearing  in  the  last  three  equations  are  defined  in  equations 
(4.1). 


4.4  PROCESSING  EQUATIONS  FOR  TRANSFORMING  ORBITAL  ELEMENTS  TO  CARTESIAN 
COORDINATES 

The  SIP  function  uses  the  inertial  Cartesian  position  and  velocity  vectors  to 
compute  inview  periods.  The  inertial  position  and  velocity  are  obtainable 
from  the  osculating  orbital  elements  generated  via  the  methods  described 
above  through  application  of  the  following  transformations: 
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4.4  PROCESSING  EQUATIONS  FOR  TRANSFORMING  ORBITAL  ELEMENTS  TO  CARTESIAN 


COORDINATES 


a(l  -  e  cos  E)  ..  1 

’ 

\yp~TF 

cos  E  / 

I  cosh  cosg  -  sinh  cos*,  sing 

-  cosh  sing  -  sinhcos*  cosg 

cosg  sinh  +  cos  h  cos^sing 

cosh  cos-i.  cosg  -  sinh  sing 

\  sioi  sing 

sin-c  cosg 
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CALL  CORDEL 


TRANSFORM  r(t'OA)- 
?(t'0A)  To  a.  e.  i. 
g.  h  AT  t'oA 

T 

SET:  At  AND  DRAG 
DECAY  RATES  TO  ZERO 

r  , 

FORM: 

_(1*  o 

i  =1 . 6 

SET:  i  =  1 


I 


CALL  BRAUER 

* 

COMPUTE  ith 
OSCULATING  ELEMENTS 

a.  e,  i,  f,  g,  h  AT 

t'O  A 

I 


FORM: 

pf'1 ;  j  =  1 .  • 

...  6 

I 


(i*1l  (il  .  (il 

-Sj  -iSj-Sj  I 

: — r~  • 

COMPUTE  (i+1)th 
MEAN  ELEMENT  SET 

a",  e",  i”,  g”,  h” 

at  t'oA 

~s~ 


FIGURE  5-1.  Process  Flow  of  the  Brouwer  Mean  Element  Converter 
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SECTION  5 

THE  BROUWER  MEAN  ELEMENT  CONVERTER  (BMC) 


5.1  FUNCTIONAL  DESCRIPTION 


The  primary  tasks  performed  by  the  BMC  are  to  receive  post  orbit  adjust 
inertial  Cartesian  position  and  velocity  vectors  at  tgA>  i,e*»  rU^)  and 


r(tQA),  and: 


(i)  transform  the  osculating  inertial  Cartesian  components  to 
osculating  Brouwer  orbital  elements; 


(ii)  iteratively  solve  for  the  Brouwer  mean  element  set  associated 
*ith  the  osculating  inertial  Cartesian  components. 


The  iterative  solution  mentioned  in  (ii)  above  utilizes  the  BRAUER  subroutine 
of  the  BOG  function.  The  process  flow  of  the  BMC  function  is  shown  in 
Figure  5-1. 


5.2  PROCESSING  EQUATIONS 

The  post-orbit  adjust  osculating  inertial  Cartesian  vectors  r(tgA)  =  (x,  y, 

z)  and  r(t qa)  =  (x,  y,  z)  are  transformed  to  osculating  Brouwer  orbital 
elements  using  the  following  relationships: 


a  = 


2 

r 


-1 


(5.1) 


4\2 

(r  •  r) 


ua 


r  r 


-  1 


(5.2) 


l  -  tan 


-1 


[  •  •  2  •  •  2  1  ^ 
|(yz  -  zy)  +  (zx  -  xz)  [ 


xy  -  yx 


(5.3) 


tad 


ITO 

vA 


w 
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5.2  PROCESSING  EQUATIONS  (con't) 


l  =  tan 


•1  I  (xy  -  yx)  y 


X  | xy  -  yx 


(5.4) 


h  =  tan" 


yz  -  zy 


xz  -  zx 


-*■  ■+ 

r  •  r 


z  r  x  r 


r|  lllid - 1  -  e2j  +  [x(zx  -  xz)  -y(yz  -  zy^l  -  e2 


g  =  tan" 


zlr  x  rlvl  -  e2  r-  * , ,r 


£  -  [x(zx  -  xz)  -  y(yz  -  zy)| [lr^-  -  1  -  e2 


(5.5) 


The  iterative  technique  used  to  find  the  mean  Brouwer  elements  from  the  as¬ 
sociated  osculating  elements  is  similar  to  that  described  by  Walter5.  However, 
a  nonsingular  orbital  element  set  devised  by  Cohen  and  Hubbard6  is  used  here 
to  provide  computational  stability  for  near  circular  orbits.  This  nonsingular 
mean  element  set,  represented  by  3y,  where 


5  Walter,  H.G.,  "Conversion  of  Osculating  Elements  into  Mean  Elements",  The 
Astronomical  Journal,  Vol.  72,  No.  8,  1967,  pp.  994  -  997. 


6  Cohen,  C.J.  and  Hubbard,  E.C.,  "A  Nonsingular  Set  of  Orbit  Elements",  The 
Astronomical  Journal,  Vol.  67,  No.  1,  1962,  pp.  10  -  15. 


? 

I 

i 
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5.2  PROCESSING  EQUATIONS  (con't) 

k  /a."\  (h"+9“+r\  .  , 

p*  cos  [y]  cos  ^ ^ - J  ,  J  =  1 

V  .  /h"  -  g“  -  l"  \  .  „ 

p  sin\T/cos\ - 2 - )  ’  J  =  2 


p^  si 


in(y)sin( 


p'4  cos 


GW 


2 

h"  -  9" 

-  r 

2 

h"  -  g" 

-  r 

2 

h"  +  g" 

+  r 

)■ 


i  =  3 


(5.6) 


e" cos  r 
-e"  sin  r 


)  •  ^  =  4 

.  J  =  5 

.  y  =  6 


is  obtained  from  the  iterative  process  executed  according  to  the  scheme 

,  /  *  1,  2,  ....  6.  (5.7) 


0  U+1)  _  &  U)  +  (£*  .  g,U)) 
Pi  J  v  i  J  ' 


The  p  in  equations  (5.6)  is  the  mean  semi-latus  rectum.  The  By  in  the  last 
equation  are  the  initial  nonsingular  osculating  elements,  and  the  3y^  are 
the  osculating  nonsingular  elements  formed  using  the  osculating  elements  ob¬ 
tained  from  BRAUER  on  the  iteration  and  generated  by  using  the  Brouwer 
mean  elements  obtained  from  By(^).  This  procedure  is  repeated  on  each  3^- 
until  the  following  convergence  conditions  are  satisfied: 


|6yU+1)  -  By^|  <  Gy 


f  y  =  i,  2,  ....  6, 


■29- 


ff-.v; 
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5.2  PROCESSING  EQUATIONS  (con't) 


where  the  t-  are  small  numbers.  When  this  condition  is  satisfied,  then 


.U+l ) 


,  j  =  1,  2,  3,  ...»  6, 


(5.9) 


*  *  %, 


and 


(!•;) 

r  i 


5  6 


(5.10) 


(5.11) 


& 


a"  =  p/( 1  -  e"2) 


(5.12) 


•V1 


i"  =  2  tan 


-l 


(J*  *  B,2)1* 


(5.13) 


3 


l"  -  tan-1  (  6s/ff  ) 

B5 


fi"  =  tan-1  (B3/-  )  +  tan-1  (BV«  )  , 


and 


Pi' 


9 "  =  2  tan"1  (BV-s- )  -  U"  +  h") 

pi 
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5.2  PROCESSING  EQUATIONS  (con't) 

For  completeness,  the  BMC  function  computes  the  Kaula  semi -major  axis  a 

i  K 

associated  with  the  post-orbit  adjust  mean  semi-major  axis  a"(t0A).  This 
is  generated  recursively,  using  the  Newton-Raphson  method,  i.e., 


*KU)  -  ^ ) 

where 

clk  +  4 5aK5  -  Y a*  +  462aK3  +  2  6  Ya^2  -  y62 

Q  ( ^  )  —  d.  —  .  ?  o 

K  K  7a*  +  206a;  -  4ya3  +  126  a";  +  46Ya^ 

K  K  K  K 


(5.17) 


(5.18) 


In  the  last  equation. 


3J  2  ( 1  -  ~2~  sin2t") 
4(1  -  e"2)3/ 2 


This  iterative  procedure  is  continued  until 


a 

ic 


U+D 


-5 

<  10 


(5.19) 


(5.20) 


(5.21) 


-31- 


NSWC  TR  82-387 


DISTRIBUTION 


Defense  Technical  Information  Center 
Cameron  Station 
Alexandria,  VA  22314 


Library  of  Congress 

ATTN:  Gift  and  Exchange  Division  (4) 

Washington,  DC  20540 

National  Aeronautics  and  Space  Administration 

Scientific  and  Technical  Library 

Code  NHS  22,  Rm.  BA39 

600  Independence  Avenue,  SW 

Washington,  DC  20546  (2) 


Defense  Mapping  Agency 
ATTN:  Mr.  Jack  Calender 
Washington,  DC  20305 


Defense  Mapping  Agency 
Hydrographic/Topographic  Center 
ATTN:  Mr.  Patrick  Fell 

Washington,  DC  20390 

Defense  Mapping  Agency  Aerospace  Center 
ATTN:  Dr.  Robert  Ballew 
St.  Louis,  MO  63118 

Naval  Electronics  Systems  Command 
Navy  Space  Project,  PME106 
Washington,  DC  20360 


Office  of  Chief  of  Naval  Operations 
Naval  Oceanography  Division  (NOP-952) 
Bldg.  1,  U.S.  Naval  Observatory 
Washington,  DC  20390 

Office  of  Naval  Operations 

Navy  Space  Systems  Division  (NOP-943) 

Washington,  DC  20350 

Naval  Research  Laboratory 
ATTN:  Mr.  A1  Bartholomew 
Washington,  DC  20375 

Naval  Oceanographic  Office 
Bay  St.  Louis,  MS  39522 


si  vx#po*v5. 


NSWC  TR  82-387 


DISTRIBUTION  (Continued) 


Office  of  Naval  Research 
Physical  Sciences  Division 

800  N.  Quincy  St.  (2) 
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Air  Force  Geophysics  Laboratory 
Hanscom  Field 

Bedford,  MA  01731  (2) 

Goddard  Space  Flight  Center 

ATTN:  Dr.  David  Smith  (1 ) 

Greenbelt,  MD  20771 

The  University  of  Texas  at  Austin 

ATTN:  Dr.  Byron  Tapley  (1 ) 

Austin,  TX  78712 

Applied  Research  Laboratory 
University  of  Texas 

ATTN:  Dr.  Arnold  Tucker  (5) 

Austin,  TX  78712 

Physical  Sciences  Laboratory 
New  Mexico  State  University 
Box  3  -  PSL 

ATTN:  Dan  Martin  (3) 

Las  Cruces,  NM  88003 

Applied  Physics  Laboratory 
Johns  Hopkins  University 
Johns  Hopkins  Road 

ATTN:  Harold  Black  (3) 

Laurel,  MD  20810 
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